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In a recent submission to this archive |lj, Suslov has claimed that our recent numerical estimate 
of the critical exponent of the Anderson transition v = 1.57 ± .02 is in error and that the available 
numerical data are consistent with a value of = 1. We contest this claim and demonstrate using 
Suslov's own scaling procedure that ^ 1. 



I. OUR PUBLISHED ESTIMATE 

In Q we published an estimate of the critical exponent 
of the Anderson transition based on a finite size scaling 
analysis of the localisation length A of electrons on quasi- 
Id bars of dimensions L x L x L^. In our published work 
L ranged from L = 4 to L = 14 and was of the order 
of 10^ in units of the lattice spacing of the tight bind- 
ing model which was studied. In more recent numerical 
studies we have extended the maximum system size for 
which data is available to L = 18. 

The critical exponent is estimated by fitting data for 
the localisation as a function of system size L and the 
magnitude of the random potential fluctuations mea- 
sured by a disorder parameter W to a finite size scaling 
form that contains one relevant scaling variable x ^^'^ 
one irrelevant scaling variable Tp, 



A = Fo (xL 



yjLyPi (xL 



(1) 

The functions Fq are F\ are expanded in Taylor series, as 
are both scaling variables. An important point to note is 
that the relevant variable is zero at the critical disorder 
Wc, so that, 

X - felT + 52t2 + . . . (2) 
while the irrelevant variable is not, so that 



V' = Co -I- ClT -f 



where 



(3) 



(4) 



In practice, when fitting our data for the box distributed 
random potential for example, we truncate the expansion 
for X 3't second order, the expansion for ip at zero order 
and the expansions for and Fi at third order. 

Using these equations to fit our numerical data we esti- 
mated V ~ 1.57±.02, where the error is a 95% confidence 
interval, for a box distributed random potential. We ob- 
tained estimates consistent with this for Gaussian and 



Lorentzian distributed random potentials. This value of 
the exponent is also consistent with a completely inde- 
pendent estimate based on the scaling of the conductance 
oi L X L X L three dimensional systems for L < 16 
and with the published estimates of other authors 



II. THE VALIDITY OF CONFIDENCE 
INTERVALS 

The fitting of a set of parameters X given a set of data 
D and a model M is often based on finding the maximum 
of the probability 



p{D\X,M) 



(5) 



This is called a maximum likelihood method since the pa- 
rameters are chosen so as to make the probability of the 
observed data a maximum. If the probability is calcu- 
lated by assuming that differences between the observed 
and predicted data are due to random measurement er- 
rors which are independently distributed according to a 
Gaussian distribution we are led to the well known x^ 
fitting procedure. 

It is crucial to note that analysis is predicated on the 
assumption that the model is correct. In particular, the 
estimate of the accuracy of the fitted parameters are all 
contingent on this assumption. To be sure, if the good- 
ness of fit is low, then one can doubt the model. However, 
a reasonable goodness of fit is not a guarantee that the 
model is correct. 

What then are the assumptions underlying our model? 
We take for granted that the Anderson transition is a con- 
tinuous phase transition which can be described by the 
renormalisation group. In addition we assume that the 
deviations from scaling in our numerical data can be de- 
scribed in practice by a single irrelevant variable. Given 
these assumptions we arrive at the accuracy quoted for 
the exponent. 

We would like to emphasise that the procedure for ar- 
riving at the accuracy of the exponent is a rational and 
rigorous procedure and the results have a precise scien- 
tific meaning and are not a matter of opinion. As with 
any mathematical procedure the validity of the results 
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depends on the correctness of the assumptions made. 
These have been clearly stated. 

It is of course always possible that if data for much 
larger system sizes were available some limitation of the 
model we use to fit the data might be exposed. If this 
turns out to be the case then our estimate of the expo- 
nent might have to be revised. Nevertheless, we strongly 
object to Suslov's characterisation of our results as "evi- 
dent disinformation." They are nothing of the sort. 

We think it worth pointing out that with the methods 
currently used for estimating the exponent the computa- 
tional time increases as L^. It seems unlikely that accu- 
rate data for substantially larger systems will be available 
any time soon. In this context we are forced to do the 
best we can with the available data by proceeding in a 
logical manner. In our opinion the libelous and ill in- 
formed comments in Rcf. contribute nothing. 

III. SPIN GLASS, WHAT SPIN GLASS? 

Suslov claims that out fit is highly non- linear and that 
maximising the probability is somehow akin to find- 
ing the potential minimum in a spin glass. Though our 
model does contain many parameters many of them are 
linear. For example, of the 12 parameters required to fit 
the data for the box distributed random potential in Ref. 
7 are linear. Hardly a spin glass ... 

IV. AN ERROR IN SUSLOV'S PAPER 

Suslov proposes what he calls "a simple procedure to 
deal with corrections to scaling." Unfortunately, his pro- 
posal is flawed because of an elementary error in Eq. (6) 
of his paper. Suslov assumes incorrectly that the irrele- 
vant scaling variable is zero atW = Wc or r = 0. In fact, 
only the relevant variable changes sign at Wc as is clear 
in the correct expansions (||) and (||) given above. If we 
look at Eq (9) of Suslov's paper we see that the quantity 
being fitted must be independent oi L a.t W — Wc- How- 
ever, the most important correction to scaling which is 
present in our numerical data is precisely a size depen- 
dence of A at Wc- Thus Suslov's method is a non- starter 
as far as modeling our data over the full range of systems 
sizes is concerned. To use Suslov's method we are forced 
to discard any data for which corrections to scaling due 
to an irrelevant variable are statistically significant. 



A = ao + Tf{L)+a2{Tf{L)f (6) 

The fitting parameters are Wc, ao, 02 and one fi = f {Li) 
for each system size present in the data. In this case this 
is a total of 8 parameters, 6 of which are non- linear. 
The results are shown in Figure |^. The critical exponent 
is then estimated by fitting the fi versus Li to a power 
law which introduces a further two fitting parameters, 
giving 10 in total. The result is = 1.53(46,60) and 
is shown in Figure ^ The numbers in brackets give the 
95% confidence interval. For comparison we also plot a 
line corresponding to = 1. We leave it to the reader 
to speculate on whether or not he or she thinks v — 1 
might be recovered for much larger systems sizes. We do 
not see any evidence for such a claim. 



VI. SUSLOV VERSUS SLEVIN- OHTSUKI 
FITTING 

When we compare the two fitting procedures the fol- 
lowing points present themselves. First, Suslov's method 
cannot describe the most important corrections to scal- 
ing which are present in the numerical data. Second, 
given that one non- linear parameter is needed for each 
system size, we end up with more, not fewer, non- linear 
parameters than the fit we used in Ref. (|]. 

VII. FINAL THOUGHTS 

In conclusion, Suslov's fitting scheme [1] does not cor- 
rectly take into account the most important corrections 
to scaling in our data. If we restrict ourselves to the 
larger systems simulated where corrections to scaling are 
negligible and apply the fitting scheme suggested in Ref. 
[1], we find v — 1.53 ± .07, consistent with our previous 
estimate in [2] but not with v = \. 

Someone once said that the only certainties in Life are 
Death and Taxes. Certainly current estimates of the crit- 
ical exponent of the Anderson transition have not reached 
this level of certainty and it is important to keep an open 
mind about how the estimates might need to be revised 
when new data becomes available. Yet we feel there is 
no need to be overly pessimistic concerning the accuracy 
of our current estimates. 



V. SUSLOV'S PROCEDURE WITH OUR DATA 



By restricting the ranges of system sizes and disorder, 
we have been able to fit some of our numerical data for 
the box distributed random potential to the following 
form based on the suggestion in Suslov's paper. 



[1] I. M. Suslov, |arXiv:cond-mat/0105325 
[2] K. Slevin & T. Ohtsuki, Phys. Rev. Lett. 82, 382 (1999). 
[3] K. Slevin, Peter Markos and T. Ohtsuki, Phys. Rev. Lett. 
86, 3594 (2001). 



2 



[4] A. MacKinnon, J. Phys. Condons. Matter 6, 2511 (1994). [7] F. Mildc, R. A. Romer and M. Schreiber. Phys. Rev. B 
[5] P. Cain, R. A. Romer and M. Schreiber. Ann. Phys. 61, 6028 (2000). 

(Leipzig) 8, SI-33 (1999). [8] F. Milde, R. A. Romer, M. Schreiber, and V. Uski. Eur. 
[6] P. Markos. J. Phys. A: Math. Gen. 33 L393 (2000). Phys. J. B 15, 685 (2000) 



3 



0.70 



0.65- 



II 

< 



0.60- 



0.55- 



0.50 




16.4 16.6 

Disorder W 



FIG. 1. A fit of our numerical data for the box distributed random potential with disorder 16 < < 17 and 10 < L < 18 
to Eq. (|). 




FIG. 2. The system size dependence of the fitting parameters f{L) from which the critical exponent is estimated. The best 
fit = 1.53 and for comparison a slope corresponding to u = 1 are shown. 
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